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Abstract. In this paper we present a rigorous asymptotic analysis for stochastic systems with 
two fast relaxation times. The mathematical model analyzed in this paper consists of a Langevin 
equation for the particle motion with time-dependent force constructed through an infinite dimen- 
sional Gaussian noise process. We study the limit as the particle relaxation time as well as the 
correlation time of the noise tend to zero and we obtain the limiting equations under appropriate 
assumptions on the Gaussian noise. We show that the limiting equation depends on the relative 
magnitude of the two fast time scales of the system. In particular, we prove that in the case where 
the two relaxation times converge to zero at the same rate there is a drift correction, in addition to 
the limiting Ito integral, which is not of Stratonovich type. If, on the other hand, the colored noise 
is smooth on the scale of particle relaxation then the drift correction is the standard Stratonovich 
correction. If the noise is rough on this scale then there is no drift correction. Strong (i.e. pathwise) 
techniques are used for the proof of the convergence theorems. 

Key words, white noise limits, Ornstcin-Uhlenbeck process, Kraichnan model, Wong-Zakai 
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1. Introduction. Many physical systems are subject to either additive or mul- 
tiplicative noise. The dynamics of such systems are quite often adequately described 
by systems of stochastic differential equations. There are various applications where 
the noise in the physical system under investigation has a non-trivial spatio-temporal 
structure and where it is not realistic to model it is as a white noise process. The 
term colored noise is used for such a noise process. 

It is a well known result that if we approximate white noise by a smooth, colored 
process, then at the limit as the correlation time of the approximation tends to zero, 
the smoothed stochastic integral converges to the Stratonovich stochastic integral 
0] 021 , [21 ch. 10]. To be precise, consider the SDE (written here in one dimension 
for simplicity) 



where b(x), f(x) are Lipschitz continuous and rj{t) is a continuous mean zero Gaussian 
process with K(rj(t)ri(s)) — ie~l t_s . Then, the results of [321 imply that, as e 
converges to 0, the solution of converges weakly to X(t) which satisfies 1 
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throughout the paper we will use the notation J f(x(s))d/3(s) (respectively f(x(t))0(t)) to 
denote the Ito stochastic integral (respectively differential) and f f(x(s)) o d/3(s) (respectively 
f(x(t)) o 0(t)) for the Stratonovich stochastic integral (respectively differential). Furthermore, we 
will refer to an Ito or Stratonovich SDE depending on how we choose to interpret the stochastic 
integral in the equation. 




(1.1) 



e 



X(t) = b(X) + -f(X)f(X) + f(x)p. 



(1.2) 
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Here (3(t) denotes a standard one-dimensional Brownian motion. The term hff is 
sometimes referred to as the Stratonovich correction. This result has been extended 
in various ways, including the case of multiple Ito integrals 19, as well as linear 
1201 |HJ ED and semilinear SPDE, an d the references therein. Moreover, the 
case of infinite dimensional noise has also been studied [H1IHIE3- I* 1 the context of the 
theory of turbulent diffusion the results of the aforementioned papers are concerned 
with convergence of rapidly decorrelating in time velocity fields to the Kraichnan 
model for passive tracers |T5] . 

The main interest of this work is to extend these results to situations where inertial 
effects arc taken into account. For motivation, consider the motion of a particle with 
relaxation time r under the influence of a force field b(x) and subject to dissipation 
and colored multiplicative noise: 

.. f(x)y(t/e 2 ) 

tx = otx) — x H . (1-3) 

e 

We are interested in analyzing the limit of 11. 3|) as both r and e tend to 0. It should 
be expected that these two limits do not commute. To see this, first let e — * while 
keeping r fixed to obtain the SDE 2 

t'x = b(x) ~ x + f(x)$. 

Taking now the limit as r -> leads to the Ito SDE |23 ch. 10] 

x = b{x) + f(x)@. 

On the other hand, if we first take the limit as r — > and then let e — > 0, we end up 
with equation (|1.2|) . 

Because of this lack of commutativity, it not clear what the limiting equation 
should be as we let e and r both tend to at the same time. This is the sometimes 
referred to as the Ito-versus-Stratonvich problem in the physics literature (26| . The 
correct form of the limiting SDE and, in particular, the presence or otherwise of a 
drift correction term in addition to the limiting Ito integral the noise induced drift- is 
of particular importance in the theory of Brownian motors |2fi| . noise induced phase 
transitions [111 1151 l2"T] and the dynamics of fronts [2Hj • 

The purpose of this paper is to investigate the Ito versus Stratonovich problem 
rigorously. Let us discuss now the main results of this paper in the one dimensional 
setting. The Id version of the model considered in this paper reads 

r e^ = b( X )-±+ f{x)v{t/e2) , (1.4a) 



fj = -ai) + VA/3, (1.4b) 

where a, A, tq are positive 0(1) parameters and 7 £ (0,oo). In this paper we 
show that three possible limiting equations result, depending on the magnitude of 
the particle relaxation time relative to that of the noise correlation time, i.e. on the 

2 It is easy to check that in this case there is no Stratonovich correction to the Ito integral, because 
of the regularity of x . 
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exponent 7. In particular, for 7 £ (0, 2) we show that the limiting equation is the Ito 
SDE 

X(t)=b(X) + ^f(X)$. (1.5) 
a 

For 7 6 (2, 00) we obtain, at the limit e — > 0, the Stratonovich SDE 

X(t) = b(X) + ^f(X) /3 
a 

= b(X) + A /W /'P0 + (1.6) 

For 7 = 2 the limiting SDE can be interpreted in neither the Ito nor the Stratonovich 
sense; we obtain: 

X(t) = b(X) + —A -f(X)f'(X) + ^f(X)P, (1.7) 

Let us make some remarks concerning (|1.7f) . If we define the stochastic integral 

* f(X(s))od/3(s) :-^77T T t f(X(s))f(X(s)) ds + ^ f f(X(s))d(3(s), 

2a z (l + T a) J a J 

(1.8) 

then this integral obeys neither the Ito nor the Statonovich (i.e. Newton-Leibnitz) 
calculus. Let us now define the stochastic integral in (|1.8|) as the limit of Ricmann 
sums 

t rr n 

f(X(s))odp(s) + (1 " M)/(*(*i-i)) 09(»j) ~ /3(*i-0), 

3=1 











with e [0, 1]. The stochastic integral l|1.8fl corresponds to the choice [1 = 2 (i+T a a) • 
As is well known, e.g. [231 ^ ^ ne Ito integral corresponds to the choice n = 0, whereas 
the Stratonovich integral to /1 = 5. Notice that letting tq in l|1.7|l vary in (0, 



00 



interpolates between these two well known integrals. 

The one dimensional model l|1.4fl was studied by Graham and Schenzle in ^2] us- 
ing formal singular perturbation analysis for the corresponding Fokker-Planck equa- 
tion in the spirit of 0]. Similar questions to the one studied in this paper were 
investigated by Givon and Kupferman in ^2] f° r SDE limits of discrete dynamical 
systems with scale separation. Specific examples were presented where the limiting 
SDE are neither of Ito nor of Stratonovich type. A formal derivation of the results 
reported in this paper together with extensive numerical simulations were presented 
in [TH]. 

In this paper we base our rigorous derivation of the limiting SDE for the infi- 
nite dimensional version of (|1.4|l -see equation below- using the pathwise techniques 
developed in [23 , following the work of Our method enables us to treat the infi- 
nite dimensionality of the noise in a rather straightforward way and, in addition, to 
prove strong convergence results. Furthermore, we are able to prove upper bounds 
on the convergence rate in L p (fl; C([0, T]; R)). The numerical results reported in [l~8*j 
indicate that the upper bounds are in fact sharp. 

This paper is organized as follows. In section[21we describe the equations that we 
will be studying and we present the convergence theorems. In section [21 we present 
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various preliminary results which are necessary for the proof of our convergence the- 
orems. In section 0] we show that the structure of the limiting equations depends 
crucially on 7. Our convergence theorems are proved in section [SJ In section we 
present two applications of the convergence theorems, with particular emphasis on 
the inertial particles problem considered in |241 1291 ISTI) . Finally, section [7| is devoted 
to some concluding remarks. 

2. Description of the Model and Statement of Main Results. In some 
of the applications of interest to us the driving colored noise is infinite dimensional. 
This arises, for instance, in Gaussian random field models of turbulence such as those 
pioneered by Kraichnan and generalizations to include noise correlation times 
[HI 171 HT H 1171 1291 Such applications are described in section H3 In this section 
we formulate the problem for infinite dimensional driving noises, and state our main 
results. 

2.1. The Model. We consider the Langevin dynamics for a particle moving in 
M. d , d > 1 under the influence of a forcing term b{x) and a rapidly decorrelating in 
time random field v(x,t): 

e l T x = b(x) + V ^ X,t ^ -x, x£R d , (2.1) 
e 

where 7 € (0, 00) and e <C 1. The field v(x,t) is a generalized Ornstein-Uhlenbeck 
process. This is a mean zero, Gaussian process which can be constructed as the 
solution of the vector valued SPDE 

dv = -Avdt + dW. (2.2) 

Here we take A : D(A) — > (L 2 (S!)) d , where fl C R d and W is a Q- Wiener process on 
H = (L 2 (n)) d . We assume that A is a strictly positive self-adjoint operator on the 
Hilbert space H and that, furthermore, it has the same eigenfunctions {fk}kLi a s Q- 

Af k = a k fk, Qfk = Xkfk- 

We now assume that there exist vectors hk € R d and positive definite selfadjoint 
operators A, Q on L 2 (Vl) such that 

fk = h k 4>k, A(j) k =a k 4>ki Q4>k — Afc^fc. 
Using this we can write 



v(x,t) = f(x)rj(t) = ^2h k (f>k{x)f]k(t), 

k=\ 

where 77(f) : £2 K is defined through the equation 

drj = -Ar] dt + dW. (2.3) 
Here, abusing notation, we have used A, Q € £(^2) with 



A = diag{afc}, Q = diag{A fc }. 
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Furthermore, W is a Q— Wiener process on £ 2 : 

oc 
k=l 

with {e-k}kLi being the standard basis in £2 and (3 k (t) mutually independent standard 
one-dimensional Brownian motions. We remark that, for each fixed x, f is a linear 
operator from £ 2 to R d : f e L(e 2 ,R d ). 

Using now the fact that f3(ct) — ^Jc(3(t) in law, we can finally write our model in 
the following form: 

e~<x = b(x) + - x, (2.4a) 



v(x,t) = f(x) V (f), (2.4b) 



di] = -^Arjdt + -dW. (2.4c) 



To simplify the notation we have set To = 1 in (|2.4a|l . In the sequel we will use the 
both notations v(x,t) and f{x)rj{t) for the random field. 

2.2. Statement of Main Results. Our goal now is to obtain the limiting 
equations of motion, as e — > 0. In order to prove our convergence theorems we 
will need to impose various conditions on the spectrum of the Wiener process, the 
eigenvalues of the operator A, the eigenf unctions {^(a;)}^! and the drift term b(x). 
The conditions that we have to impose are more severe for 7 > 2, since in this 
parameter regime we will need more integrations by parts in order to obtain the 
limiting equations. 

We will use the notation || • || to denote the Euclidean norm in R d . Subscripts 
with commas will be used to denote partial differentiation. 

As has already been mentioned, we take A to be a self-adjoint, positive operator 
on L 2 (tt). We assume that the eigenvalues {a k } k x =1 of A satisfy 

• • • > afc+i > a k > u > 0, V tt^ < (2-5) 

2«fe 

fc=i 

The eigenfunctions of A are normalized so that their L 2 {Q) norm is set to 1: ||</>fc||z, 2 (n) = 
1. Moreover, for 7 S (0,2) we assume that there exist constants C > 0, a, /? such 
that 

r M*)£C 2 (n), k = i,2,..., 

I ||^fcWlk~(n) < Ca^, \\D<pk(x)\\L<»(n) < C a p k . 

The conditions for 7 > 2 are more severe. We assume that there exist constants 
C > 0, a, /?, 7, S such that 



fc = l,2 

«g, \\D 

\D 2 cp k {x)\\ L ~ (n) < Cal \\D 3 Mx)\\L~(n) < C a{ 



k(aOIU-(n) < Cat, ||X>^fc(x)||L-(n) < Cof, (2.7) 
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Remark 2.1. At this level of generality and, in particular, since we do not make 
any specific assumptions on the operator A, we do not have any detailed information 
on the L°° norm of the eigenf unctions {<pk\kLi and their derivatives. Much is known 
when A is a uniformly elliptic operator, see e.g. L'Ml Ch. 5], U\\14Jj and the references 
therein. In particular, the results from \14\ imply that, when A is a uniformly elliptic 
operator with smooth coefficients and Dirichlet or Neumann boundary conditions on 
some bounded domain fl C K rf with smooth boundary, then the following estimate 
holds: 

\\D n <f> k \\ Loam <Ca^-,n = 0,l,..- (2.8) 

We will assume that the drift b(x) is Lipschitz continuous: 

\\b(x)-b(y)\\<C\\x-y\\, x,yeR d . (2.9) 

Moreover, we will assume that there exist constants C, r such that 

\\h k \\ <C|a fe | r , fc = l,2,...,. (2.10) 

Now we are ready to present the conditions that we have to impose on the spectrum 
of the Wiener process. First, we need to ensure the existence and uniqueness of the 
equations of motion l|2.4a|l . To this end, we assume that the velocity field is sufficiently 
regular: 3 

v{x,t) e (C{R + ,c 1 (n))) d . (2.11) 

Assumption l|2.11|L together with assumption l|2.9|l ensure that there exist almost 
surely a unique solution of the equations of motion (|2.4a|) . when the initial condi- 
tions for l|2.4c[) are distributed according to the invariant measure of this process. 
Furthermore, for 7 < 2 we have to assume conditions of the form: 

00 

]T VA^ir^^ <oo, (2.12a) 

k=\ 
00 

J2V^4 +fJ -i- p) <oo, (2.12b) 
fe=i 

3 A simple variant of gSJ Thm. 5.20] yields that v{x,t) G (C(R+ , C 1 (T2))) d provided that there 
exists 8 ( 6 (0, 1) such that 

00 

E, 2(r+a)-l-C , 

fc = l 

00 

E, 2(r+C7)-l ^ 

fc = l 

However, these conditions are not optimal and so we simply assume 12.111 . 
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The specific value of the exponent p will be given when stating our convergence 
theorems. For 7 > 2, in addition to the (|2.12|) we further assume: 



^ /Ajtafc 2 < 00 (2.13a) 
fe=i 



E 

fc=i 



A fc aj. < 00 



(2.13b) 



k=l 



\, 2r+/3+ 7 - 



< OO 



(2.13c) 



Remark 2.2. Consider the case where A is a uniformly elliptic operator. From 
(|2.8(l it is easy to see that conditions 12.12J1 and 12.13(1 become 

00 

E ^r d/4 - 5/8 - p) < 00, 
fc=i 

and 

00 

E, 2r+d/2-17/8 . 
A fe a fe <oo, 

fc=l 

respectively. We fix now an integer p > 1. We assume that the initial conditions are 
random variables, independent of the er-algebra generated by W(t), with 

E||a: || 2p < oo, E||yo|| 2p <oo, (2.14) 

for 7 e (0, 2) and 

E||:r || 2p < 00, E||iyo|| 4p < 00, (2.15) 

for 7 S [2, 00). 

Now we state the convergence theorems. We start with 7 € (0, 2). 
Theorem 2.1. Let x(t) be the solution of equation (|2.4a|l for 7 e (0,2). As- 
sume that conditions <j!TfJj) . (|23|) . ljTTU|l . (|2~TT1) . u«#i p = | and ifTTljl 

ZioZd. Assume further that the initial conditions for l|2.4c|) are stationary. Then x(t) 
converges, as e — > 0, to X(t) which satisfies the following equation: 

X(t)=x + [ b(X(s))ds+ [ fiXis^A- 1 dW{s), (2.16) 
Jo Jo 

the convergence being in L 2p (Q, C([0, T],R)): 

E ( sup \\X(t) - x(t)f p ) < C (e<v + e( 2 -^P- a ) , (2.17) 

where a > is arbitrarily small. The constant C depends on the moments of the 
initial conditions, the spectrum of the Wiener process, the operator A, the exponent 
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p, the maximum time T and a. In order to present the convergence theorems for 
the case 7 > 2 we need to introduce some notation. We denote by 0, 9 : £2 — ► ^2 the 
diagonal operators defined by 



6 = diag^^ e = diag<j ^ 2 ,/ i — \. (2.18) 



We will use the notation V • A to denote the divergence of a matrix A, i.e. {V • A}i = 
The next theorem covers the case 7 S (2, 00). 

Theorem 2.2. Letx{i) be the solution of equation <|2.4a[l /or 7 S (2, 00). Assume 
that conditions l|23)l. (|2~7|l . (|2~U|) . lj2~TU|l . (I?!TT|) . I|2.12a|) imi/i p = \, (|2.12b|) im'tt 
/? = 0, (|2.13() and (I2.15|) /10W. Assume further that the initial conditions for (|2.4c|l 

are stationary. Then x(t) converges, as e — > 0, to X(t) which satisfies the following 
equation: 



t rt 



X(t) = x Q + / b(X(s))ds+ / V • (f(X( S ))ef T (X(s))) ds 
Jo Jo 

f(X(s))OV ■ f T (X(s)) ds + f fiX^A- 1 dW(s), (2.19) 

Jo 

the convergence being in L 2p (fl, C([0, T],R)): 







E 



( sup \\X(t) - x{t)\\ 2 A < C (e 2p - a + e M~f-i)-A t (2.20) 

\0<t<T J V ' 



where a > is arbitrarily small. The constant C depends on the moments of the 
initial conditions, the spectrum of the Wiener process, the operator A , the exponent 
p, the maximum time T and a. Finally, the case 7 = 2 is covered by the following 
theorem. 

Theorem 2.3. Let x(t) be the solution of equation (|2.4a|l for 7 = 2. Assume 
that conditions (ESJ, $Fft, l(2~9")). (jlHUj) . (|2~TT|) . (|2.12al) wtfi p = ~, (|2.12b|) wit/i 
/? = 0, (|2.13|) and (I2.15|) /10W. Assume further that the initial conditions for <|2.4c[l 

are stationary. Then x(t) converges, as e — > 0, to X(t) which satisfies the following 
equation: 

X(t) = *o + £ b(X(s)) ds + j\- (/(X( S ))6/ T (X( S ))) ds 

f(X(s))eW-f T (X(s))ds+f fiX^A-UWis). (2.21) 

Jo 

the convergence being in L 2p (fl, C([0, T],R)): 

e( sup ||X(t)-s(t)|| 2p ) <Ce 2p -°, (2.22) 

\0<t<T / 

where a > is arbitrarily small. The constant C depends on the moments of the 
initial conditions, the spectrum of the Wiener process, the operator A , the exponent 
p, the maximum time T and a. 

Remark 2.3. The second and third integrals in ()2.19|) give the d-dimensional ana- 
logue of the Stratonovich correction ^f(X)f'(X) in (|l.(j|) . when the system is driven 
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by an infinite dimensional noise process. Similarly, the second and third integrals in 
(I2.21|) correspond to the drift correction in equation (|1.7fl . 

Remark 2.4. The assumptions of the convergence theorems ensure Lipschitz con- 
tinuity and linear growth of all terms that appear in the limiting equations and, hence, 
existence and uniqueness of solutions. 

Remark 2.5. Throughout the paper we have set To = 1, in order to simplify the 
notation. Of course, the above convergence theorems hold true for arbitrary To > 0. 
In this case, the matrix widehatO defined in i|2.18f) has to be modified: 

6 = diag \ n 9 ./ j 1 . (2.23) 

Notice that we can formally retrieve the limiting equation for 7 < 2 and 7 > 2 by 
sending tq in (|2.23|l to 00 and 0, respectively. 

Remark 2.6. Problem (I2.4|l for 7 = was considered in It was shown there 
that, under appropriate conditions on the spectrum of the Wiener process and the 
operator A, the particle position x{i) converges pathwise to the solution X of a second 
order SDE which we formally write as 

X = b{X) -X + f{X)A~ 1 W. 

It was proved in that the convergence rate is ofO(e 2 ~ a ), where a > is arbitrarily 
small. It is natural, therefore, that the convergence rate in Theorem \2.1\ degenerates 
as 7 tends to either § or 2, since the limiting equation is different in both cases. 

2.3. Remarks on the Convergence Theorems. We present now a few com- 
ments on the convergence theorems. First, we note that the smoothness assumptions 
on the eigenfunctions {4>k}kLi are more severe for 7 S [2, 00) than for 7 £ (0, 2). This 
is because, in order to prove our convergence theorems for 7 > 2, we need additional 
integrations by parts, using the ltd formula. As a result, we need to assume that 
more moments of the particle velocity at time t = exist when 7 > 2. Notice also 
that the convergence to the limiting equations becomes arbitrarily slow as 7 — ► and 
7 — * 2~ in Theorem l2.ll as well as 7 — * 2 + in Theorem l2.2l This is also not surprising 
since the form of the coefficients in the limiting equation is discontinuous at 7 = 2. 
The extensive numerical experiments reported in 18 indicate that the convergence 
rates of our Theorems are sharp. On the other hand, the conditions that we have to 
impose on the spectrum of the Wiener process, conditions l|2.12|l and l|2.13|l . are not 
sharp and are not-in general-independent from one another. In order to optimize 
these conditions one needs more detailed information on the specific problem under 
investigation, in particular on the properties of the eigenfunctions of the operator A. 
Consider for example the case where A is a uniformly elliptic operator, Remarks 12. II 
andl2~21 

Let us now try give an intuitive explanation of our results. First, for 7 < 2 
the particle relaxation time — which is of C(e 7 ) — is large compared to the relaxation 
time of the noise-which is of C(e 2 )- and consequently the particles experience a rough 
noise with practically zero correlation time. This means that for 7 < 2 the OU process 
is not viewed from the point of view of the particle as a smooth approximation to 
white noise and, as a result, the stochastic integral in the limiting equation has to be 
interpreted in the ltd sense. On the other hand, when 7 > 2, the particle relaxation 
time is small compared to that of the noise. Consequently, in this parameter regime 
the rescaled OU process is indeed a smooth Gaussian approximation to white noise and 
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the stochastic integral in the limiting SDE should be interpreted in the Stratonovich 
sense as in equation <|2.19[1 . in agreement with standard theorems [21 sec. 10.3]. The 
case 7 — ► oo leads to the case of tracer particles whose relaxation time is zero and 
covered precisely by these standard theorems. 

For the case 7 = 2 the particle relaxation time is comparable in magnitude to 
the noise correlation time and a resonance mechanism prevails which results in the 
limiting stochastic integral being neither that of Ito nor that of Stratonovich. In 
this case the drift correction to the Ito stochastic integral depends on the detailed 
properties of the OU process, in particular its covariance. 

It is well known that for second order stochastic differential equations the Ito and 
Stratonovich interpretations of the stochastic integral coincide. For certain Gaussian 
fields vix, t) this also happens for the limiting equations given in our convergence 
theorems: the Stratonovich correction, as well as its modified version from Theorem 
12.31 will, in some situations, be identically zero due to the specific properties of v(x, t). 
In this case the limiting equations are the same for all values of 7. This situation occurs 
for example in the inertial particles problem which is discussed in section due to 
the fact that the fluid velocity is assumed to be homogeneous and incompressible. 

Let us now outline the method that we will use in order to prove the results of 
this paper. The first step is to use the variation of constants formula to write the 
particle velocity y(t) := x(t) and particle position x(t) as follows: 



respectively. The next step is to use equations (|2.24|l and (|2.25|l in order to obtain 
sharp estimates on the moments of the particle velocity. The basic strategy will be to 
derive first estimates valid for 7 £ (0, 00) and then use them in order to obtain sharper 
estimates valid for 7 g (0, 2). We emphasize that sharper estimates for 7 S (0, 2) are 
necessary for the proofs of the convergence theorems. Now, with the estimates for 
the moments at hand we prove that the last two integrals on the right hand side of 
equation fFZty are small in L 2 *>(ft; C([0, T], M)) for all values of 7 > 0. 

Then we study the term which induces noise in l|2.25[l . in the limit e —> 0, namely 
the first integral on the right hand side of this equation. We refer to this as I(t). 
We use the Ito formula, together with the estimates on the moments of the particle 
velocity, to show that I(t) consists of an 0(1) term plus higher order corrections. The 
leading order term in I(t) is different for 7 < 2, 7 > 2 and 7 = 2: this is the term 
which is responsible for the difference in the limiting equations for different 7. Finally, 
the proof of the convergence theorems is completed by an application of Gronwall's 
lemma. 

Throughout the paper we will make extensive use of estimates on the infinite 
dimensional OU process v(x, t) as well as the stochastic convolution 




(2.24) 



and 




(2.25) 
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Lemmas 13 . 21 and 13.41 respectively. The proof of the first of the above lemmas is based 
on Borell's inequality from the theory of Gaussian processes PQ, while the proof of 
the second uses the factorization method |25|. 

We remark that, unlike the methods used in the proofs of the convergence the- 
orems in |24| . the proof in this paper relies on the presence of the friction term —x 
in the equations of motion l|2.1|l . The linear friction term enables us to obtain repre- 
sentations H2.24fl and (|2.25() for the particle velocity and position, respectively, which 
are necessary for analyzing the dependence of various moment bounds on e. 

3. Preliminary Results. 

3.1. The Integral Formulation. The first step is to obtain an integral equation 
for x(t) that will be more convenient for our analysis. 

Lemma 3.1. Consider the equations of motion (I2.4a|) . Then the particle position 
satisfies the following integral equation: 

, x -v /, _ ' > f /„ v(x(s),s) , 

x(t) =x + e 7 y (l ~e 7r )+ 1 ~ ds 



+ (l - e^) b(x(s))ds. (3.1) 



Proof. We start by solving the equation for y(t) = x(t) using the variation of 
constants formula: 



2/ ( f )= yoe -iV +e -7 / e V ?^gWli2 ds + e -7 / e ^b(x(s))ds. (3.2) 



Another integration will give us an integral equation for x(t) which involves a double 
integral: 



x(t) = xq + / j/o e ^ ds -\ / (/ e' 71 ' H(x(s), s) ds | dl 

Jo e 7 Jo \Jo J 

= x Q + eiy (l-e-^)+I(t), (3.3) 



where 



v{x(s),s) 
H(s) = h b{x(s)). 



We can reduce I(t) to a single integral as follows: First we define the following func- 
tion: 



F(£) = [ e^H{s)ds. 
Jo 
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Now we perform an integration by parts: 
1 



\J0 \ 

t 

e-^F(£)d£ 

o 



\e=o 



o 



-tt f e ^H(s)ds+ [ H(s)ds 
Jo Jo 

(l-e&) ds + J* (l - e^) b(x(s)) ds. (3.4) 

Substituting (|3~3|) into ij^t we obtain JXTJ- C 

Throughout the paper we will use the following notation: 

h{t) = e 7 y (l - e—&), (3.5a) 



Ja(t)= rm^m^ (3 . 6b) 



Ja(t) = - [\^lMf)m ds , (3.5c) 



h(t) = -[ e^b(x(s))ds. (3.5d) 
Jo 

Using this notation the particle position x(t) can be written in the form 



4 -« 
i=l 



4 

i(t) = x + V 7i(t) + / b(x(s)) ds. (3.6) 



We clearly have: 



E ( sup ||/i(t)|| 2p ) < Ce 27p . (3.7) 



Now we want to study terms hit), hit) and hit). As explained in section l2~3l we 
want to show that I 3 (t) and I 4 (i) are o(l) in i 2p (fi, C([0,T],M)) for every 7 g (0,oo) 
and then show that the behavior of the term hit) as e — > depends on 7. In order 
to obtain the necessary bounds we will need sharp estimates on the moments of the 
particle velocity. We will obtain these estimates in section 1X51 Before doing this, we 
need some estimates on the velocity field v(x,t). 
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3.2. Estimates on the Colored Noise. In this subsection we present two 
results which will be used in the proofs of the convergence theorems. We start with 
an estimate on the infinite dimensional OU process. 

Lemma 3.2. Assume that conditions (|2.7|) . I|2.1U[1 and ()2.12afl with p e K are 
satisfied. Then the following estimate holds: 

E( sup \\A- p v(x{t),t)\\ 2p j < Ce~ a , (3.8) 

\0<t<T J 

where a > is arbitrarily small. 

Proof. Let rj(t) — {rjk(t)} k % 1 : t 2 — * K be the solution of (|2.3|l with stationary 
initial conditions. The kth component r/ k (t) solves the equation 

dn k = -a k n k dt + y/x^dWk- 
A simple variant of Theorem A.l from [21] yields: 

E bf,"*" 2 *) s fe) s '( i+ K^))- (3 - 9) 

Let Y(x(t), t) := A~ p v(x(t), t). We first consider the case p = 1. We have: 



E sup ||r(x(t),<)f = 

\0<t<T 



< 



sup 

0<t<T 



sup 

0<t<T 



^2h k cl) k (x(t))a k p r lk (t) 
fc=l 

00 

^ • hi(j) k {x{t))(j)i{x{t))al P aJ p n k {t)irii(t) 

k,t=l 

s( sup (£||M 2 l^(sW)lV 2p M*)l ; 

-E[ sup (^^||^||||^||0 fc (a;(t))^(a;W)afc P a7 P ?7 fe (t)r7 £ (t) 



0<t<T 



, fe=l ^fe 



, , |L|2(r+a) , - , 



k=l 



0<t<T 



\k\ 



r+a 



\k=l k 



SUp \Tj k 
0<t<T 



Kk=l a k 



Xk\k\ r+a 



\k=i a, 



on account of condition (|2.12a|) . We can proceed in the same way for p > 1, by 
breaking the sums into various parts, until we have sums that involve independent 
OU processes. Condition (|2.12a|l ensures the summability of all the sums that appears. 
The lemma is proved. □ 



Using the above lemma we can easily obtain the following estimate. 
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Lemma 3.3. Assume that the conditions of Lemma V) '.ifl hold. Define 

I(t) = / e^-A- p v(x(s),s)ds. 
Jo 

Then the following estimate holds: 

sup ||/(t)|| 2p ) < Ce 2p ^ a , 

0<t<T ) 

where a > is arbitrarily small. 
Proof. We have, for t e [0, T]: 



(3.10) 



(3.11) 



\I(t)\\ 2p = / e^A-Pv(x(s),s)ds 



■2p 



< sup \\A- p v{x(s),s)\\ 2p 

0<s<T 



e"^ ds 



o 



<e 2p '< sup \\A' p v{x{s),s)\\ 2p . 

0<t<T 

Lemma f3 . 21 now yields estimate (|3.11|l . □ 

Remark 3.1. The techniques used in the proof of Lemma \H.°A enable us to conclude 
that we can bound uniformly all moments of the field v{x,i): 

¥,\\A- p v{x,t)\\ 2p < C, 

provided that the assumptions of the Lemma are satisfied. Furthermore, the method 
of proof of Lemma \S.'A gives: 

E||/(t)|| 2p < Ce 2p ~<, 

where I(t) is defined in l|3.10[l . 

Remark 3.2. Lemma \S. 'A with p — provides us with estimates for Is(t) and Ia(£) 
which we will use for 7 € [2, 00): 



E 



( sup ||/ 3 (i)|| 2p ) <Ce 2p ^-\ 7 E [2,oo), 

\0<t<T J 



and 



E( sup \\h{t)\\ 2p )<Ce- 2p -°, 7 e[2,oo), 

vO<t<T 



(3.12a) 



(3.12b) 



where er > is arbitrarily small. This estimate is not sharp enough when 7 6 (0, 2) 
and we need to improve it. This will be accomplished in Corollaru \3.1(A 

Remark 3.3. Assume that the moments of the particle velocity y{t) satisfy 



E ( sup 

VO<t<T 



l|yWII 2p ) <Ce< p , 



for some (gl. A repeated use of Holder's inequality, together with the Gaussianity 
of the process rj(t) as in the proof of Lemma 4-3 in WM enables us to prove that 



E( sup \\y(t)\\ 2p \\ V (t)\\'P) <e 

\0<t<T 



(3.13) 
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for every n, assuming that Tr(Q) < oo, with a > ; arbitrarily small. In the sequel 
we will have the occasion to use estimate i|3.13[l and variants of it repeatedly. We 
proceed now with an estimate on a stochastic integral. 
Lemma 3.4. Consider the stochastic integral 

I(t)= / e^f{x(s))A- l dW{s). 
Jo 

Assume that conditions (|2.6|) . (I2.1U|) and (|2.12a[l with p — \ hold. Then we have the 
following estimate: 



El sup ||/(i)H 2p <Ce pl -° 

K 0<t<T 



(3.14) 



where a > is arbitrarily small. 

Proof. We fix a € (0, h) and use the factorization method from sec. 5.3] to 
obtain 

1(f) := f e^f(x{s))A- 1 dW{s) 
Jo 

sin(7ra) ''' 



(f- s) a - i Y{s)ds 



Y(s) = / (s - a)- a f{x{a))A- 1 dW(a). 



where 



We choose m > ^- and use Holder inequality to obtain: 



I(t)\r<C e^(t-s) 



ds 



\\Y(s)\\ 2m ds. 



A change of variables now yields: 



J(t) 



o 



e^r(t-sY 



ds 



2m- 1\ 2 



-(o+2(m-l)) 



2m 



In the above estimate we used the fact that, since m > tt-, we have e z z ! ™- 1 ^ a ^ € 
i 1 (R+). Consequently we have: 

E^sup ||/(i)|| 2m> ) < CeT^-^E / ||y(s)|| 2m cfo. 

\0<t<T / Jo 

To proceed further, we use Lemma 7.2] to deduce that there exists a constant 
C m > depending only on m such that 



sup E||F(s)|| 2m < C m E 

0<s<T 



-(s-a)- 2a \\f(x(a))A-Xoda 
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with 

Wfixia^A-^ := Tr [(f{x{a))A^) Q (/(^a))^ 1 )* 

= ^ Xk\M^))\ 2 \\hk\\ 2 
fc=i fe 



< C^Afclafel 2 ^''- 1 ) <oo, (3.15) 



k=l 



on account of condition (|2.12a|) with p = h . Now we can apply the same change of 
variables that we used in the estimate for J(t) to obtain: 

sup E||F(s)|| 2m < CE ( / e 2£ ^ {s - a)- 2a da 

0<s<T \Jo 

<C [t l{1 ~ 2a) e~ z z~ 2a dz 

< Q e 7m(l-2a) 

From the above estimates we conclude: 



E ( sup \\I(t)\\ 2m ) < Ce^™- 1 '. 

\0<t<T J 



Now estimate l|3.14ll follows for p > ^ upon taking p — m. For p < j- we apply 
Holder inequality to obtain: 

fif SUp \\I(t)\\ 2P ) < M BU P \\I{t)\\ 2m )Y <C el*-*, 
\0<t<T J \ \0<t<T J J 

which completes the proof of the lemma, since m can be chosen to be arbitrarily large. 
□ 

3.3. Bounds on the Moments of y(t). In this subsection we will obtain 
bounds on the moments of y(t) that we will need for the convergence theorem. In 
order to obtain estimates on the moments of the particle velocity we first need to 
obtain a crude estimate on the moments of the particle position. This estimate will 
be improved later. 

Lemma 3.5. Let x(t) satisfy equation \2.4a\ l- Assume that (|2.9[) and <|2.14ll as 

well as the conditions of Lemma with p — are satisfied. Then the following 
estimate holds: 

E ( sup \\x{t)\\ 2 A < Cr 2p - ff , (3.16) 

\0<t<T / 

where a > is arbitrarily small. 

Proof. The particle position is given by 1)3. ljl , which can be written in the form 

x(t) = x + h(t) + I 2 (t) + I 3 (t) + j (l - e^*) b(x(s)) ds. 
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The Lipschitz continuity, assumption l]2.9p. of b(x) implies that there exists C > 
such that 

\\b(x)\\<C(l + \\x\\). (3.17) 
We use this, together with estimates l|3.7[l and l|3.12[l as well as Lemma f3. 21 to obtain 



E 



( sup \\x{t)\\ 2p ) <E\\x \\ 2p + Ce 2 ~<P + Ce- 2p - a + Ce^-V 

\0<t<T J 

+ ( E f sup \\x(t)\\ 2pS ] dt 

Jo \0<t<T J 



-1 -CT 



+c 

<Ce- 2p - a + [ T e( sup \\x{t)\\ 2p ) dt. 

JO \0<t<T ) 

Estimate (|3~To|) now follows from Gronwall's lemma. □ 

We have already mentioned that the bounds on the moments of the particle 
velocity will be different for 7 < 2 and 7 > 2. We start with the regime 7 € [2, 00). 

Lemma 3.6. Let x(t) satisfy equation j2.4a)) an let y(t) — x(t). Assume that 
assumptions l|2. 9J) and H2.14|l as well as the conditions of Lern,ma \H . 6 A with p = are 
satisfied. Then the following estimate holds: 

e( sup \\y(t)\\ 2p ) <Ce- 2p -\ (3.18) 

\0<t<T / 

where a > is arbitrarily small. 

Proof. The particle velocity is given by l|3.2|) . which can be written in the form 

y{t)=yoe-^ -e^h{t)-e^h(t). 
Lemma f3. 51 and estimate (|3.17|) give 

E ( sup ||/ 4 W|| 2p ) < Ce 2p ^- 1) - a . 

\0<t<T ) 

We use the above estimate, together with (|3.12a|) to obtain 



E ( sup 

V0<t<T 



2 p 



\\y(t)\\ 2p ) <Cx + C 2 e- 2 ^El sup \\h{t)\\ 2p \+C 3 e-^ p E( sup \\h(t)\ 

J \0<t<T J \0<t<T 

< C(e- 2p - a + 1), 
from which the estimate follows. □ 

Estimate (|3.18|l will be sufficient for our purposes for 7 £ [2, 00). However, it is 
not sharp enough for 7 £ (0, 2). In order to prove the convergence theorem for values 
of 7 in this parameter regime we need to improve the estimate for the 2pth moments 
of the particle position and particle velocity. For these two estimates we need some 
preliminary estimates which will be also used in the proof of the convergence theorem. 
We start with the following lemma. 
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Lemma 3.7. Assume that conditions l|2.5f) . (|2.6|l . I|2.1L)[I . together with (|2.12b() 
with p £ R hold. Fix x, y <E R d and define df(x)y e L(€ 2 ,K d ) by 



{(df)y} 1 = J2 h kVV(b, 



fc7fc- 



k=l 



Let 



and 



I(t)= / e^df(x(s))y(s)A-"T ] (s)ds 



I(t)= / df(x(s))y(s)A-Pr)(s)ds. 
Jo 

Then the following estimates hold: 



E( sup 

\0<t<T 



\\I{t)\\ 2p ^j < Ce<- 2p ~ 1 ^- 2p - 



(3.19) 



(3.20) 



E ( sup 

\0<t<T 



1(f) 



2;» 



where a > is arbitrarily small. 

Proof. An application of Holder inequality yields: 

2P-1 ,t 



(3.21) 



\\I{t)\\ 2p < U ei^^'dsj J \\df(x( S ))y( S )A-P V (s)\\ 2p d S 
< e (2p-l) 7 f \\df(x(s))y(s)A-Pn(s)\\ 2 Pds. 



A simple variant of l|3.13|l . together with calculations similar to those used in the 
proof of Lemma 13.21 gives 

e( sup \\I{t)\\ 2p ) <e {2p - 1 ^ [ T E\\df{x(s))y{s)A' p r 1 (s)\\ 2p ds 

\0<t<T J Jo 



< Ce ( - 2p ~ lS> ' 1 



E 



oo d 

E hA x ( s ))yk(s)aJ p r tj (s) 
j=l k=i 



ds 



< 



( / \ 2p 

Ce^^J^ iM j ||y(,)| 



2p 



< c<P' p ~ 1 ^ 1 ~ 2p ~ a 



assuming that condition (|2.12b|) holds. This proves l|3.20|l . The proof of l|3.21|l is 

almost identical and is omitted. □ 
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We introduce some notation that we will use repeatedly throughout the rest of 
the paper. We set: 



A"' i -v{x{t),t) - e~^v{x Ol 0) , J 2 (t) := / e ^ df(x(s))y(s)A" 1 r 1 (s) ds 



Ji(i) 



J 3 (t):= / f{x{s))A~ 1 dW(s), J 4 (t) := / A- 1 v{x{s),s)ds 



Lemmas 13.21 |3~3I3.4I and 13.71 and assuming that the initial conditions for v(x,t) are 
stationary and that conditions (|2.12a|) with p = \ and i|2.12b(l with p — 1 are satisfied, 
provide us with the following bounds: 



E ( sup ||Ji(i)|| 2p ) < Ce 

\0<t<T 



El sup \\J 3 (t)\\ 2p ) <Ce p ^ a , El sup ||J 4 (i)ir ) < Ce 



( sup || J 2 (i)|| 2p ) < ce^- 1 ^-^-" ; 

\0<t<T J 



0<t<T 



0<t<T 



2p7-cr 



(3.22) 



Lemma 3.8. Let x(t) be given by 1)3. 1|1 and let y(t) — x(t) be given by 1)3. 2J1 . 
Suppose that conditions l|2.9|l and (|2.12|l with p — 1 /ioW. TTien /or 7 G (0, 2) i/ie 
following estimate holds: 



E( sup ||x(i)|| 2p ) <C. 

vO<t<T 



(3.23) 



Proof. Consider first the term I-z(t) defined in (|3.5b|) . We integrate by parts to 
obtain: 



where 



and 



h(t)= / f(x(s))A- 1 dW{s)^eJ 1 {t) + eJ 2 {t), 
Jo 

Ji(*) := [A-'vixit)^) - A^v(x{0),0)] , 



(3.24) 



J 2 (t)= I df(x(s))y(s)A- 1 r,(s)ds. 
Lemma \'S . 21 with p = 1 gives 



sup 

0<t<T 



eJi(t) 



2p 



< Ce 2p - a . 



Moreover, Lemma \'S . 71 with p = 1, in particular estimate (|3.21|) . yields 

2pN 



E ( sup 

\0<t<T 



eJ 2 {t) 



< C. 



Furthermore, the Burkholder-Davis-Gundy inequality, assuming that condition H2.12a|) 
with p = \ holds, yields 



E sup 

\ 0<t<T 



f(x{s))A- x dW{s) 



2/A 



< c. 
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We put the above estimates together to conclude that 

2" 



E sup \\I 2 (t)\r) <C. 

\0<t<T J 

We use this estimate, together with ()3.12al) and the assumptions Ij2.14|l in equation 
(|3.1|1 to obtain: 



E 



( sup \\x(t)\\ 2p ) <CE\\x \\ 2p + CE\\y \\ 2p + Ce 2p{ '*- 1) - a + C [ E f sup \\x{t)\\ z ) dt 

\0<t<T J Jo \0<t<T 

<C+ [ E f sup ||x(t)|| 2 ) dt. 

Jo \0<t<T J 

Estimate (I3.23|) now follows from Gronwall's lemma. □ 

Now we are ready to obtain a sharp bound on the moments of the particle velocity 
for 7 G (0,2). 

Lemma 3.9. Let x(t) be the solution of i|2.4a[l and let and y(t) = x(t) . Assume 
that E||yo|| 2p < °° and that conditions (|2.6[) . l|2.9l) . (|2.10|l and (I2.12al) with p= \ and 
(|2.12bl) with p ~ 1 are satisfied. Then the following estimate holds: 

E ( sup \\y(t)\\ 2p ) < Ce^ p -\ 7 € (0,2). (3.25) 



Proof. Lemma \'3. 81 and estimate (|3.17|l give 

E ( sup ||/ 4 WH 2P ) < Ce 2p7 - 

\0<t<T / 

We combine this estimate with (|3.2(l to deduce: 



(3.26) 



E ( sup 

V0<t<T 



\\y(t)\\ 2p ) < Ci + C 2 e- 2p ->E ( sup ||7 8 (*)|| 

/ \0<t<T 



2p 



We need to get a sharper estimate on J3 (t) than (|3.12a() . For this we need to integrate 
by parts. We apply Ito formula to the function 



G ik (s,x,r)) = e^fi k {x(s)) 



Vk(s) 



to obtain, after some algebra: 

I 3 (t) = e \A~ x v{x{t),t) - e~^A~ x v(x o ,0) 
e~^~ df(x(s))y(s)A~ 1 T](s) ds 

fixis^A- 1 dW(s) 



+e 



1-7 



e~ rr A 1 v(x(s),s)ds 



--: cJi(t) - eJ 2 (t) - J 3 {t) + e 1 -V 4 (t). 



(3.27) 
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Consequently, on account of estimates l|3.22|l : 

E ( sup \\I 3 {t)\\ 2 A < C (e 2p - a + e^-^T-^ + t^-*) . 

\0<t<T ) ^ ' 

Thus: 

E ( sup ||y(i)|| 2p ) < C (>p(i-7)- CT + e -i-° + e -P7-<A 5 

\0<i<T / ^ ' 

from which estimate <|3.25[) follows upon noticing that, for 7 € (0, 2), 2p(l — 7) > — P7 
as well as that p > 1 4 . The proof of the lemma is now complete. □ 

From the above lemma we can obtain sharper bound on 13(f) and J2(t) defined 
in fTTOIl for 7 € (0,2): 

Corollary 3.10. Lei x(i) 6e £/ie solution of (I2.4a|) . Suppose that the conditions 
of Lemma hold and that 7 € (0,2). Then 13(f) and J2(t) satisfy the following 
estimates 

e( sup \\I 3 (t)\\ 2p ) <Ce™-°, 7 €(0,2), (3.28) 

\0<t<T J 



and 



J 2 (t) 2P ) <Cr^, 7 G (0,2), (3.29) 



E sup 

\0<t<T 

respectively, where a > is arbitrarily small. 

Proof. Consider the term Jz(t) defined in equation (|3.27|l . Lemmas 13.71 and 13.91 
imply that for 7 € (0, 2) we have: 

e( sup ||J 2 (t)|| 2p ) < C (e^- 1 ^-* + e 2(p-D7-^ 

< Ce^- 1 ^-' 7 , (3.30) 
since p > 1. We use now (I3.27|) and the above estimate, together with (|3.22|) to obtain: 

E( sup ||/ 3 (<)|| 2p ) < C fe 2p - CT + e 2p-+(f-l)7-^ + e Pi-A 

\0<t<T J \ ' 

where we have used the facts that p > 1, 7 < 2. Calculations similar to the ones used 
in the proof of Lemma |3~7I together with estimate l|3.25|l yield estimate <|3.29[1 . □ 

4. Drift Corrections To the Ito Integral. Estimates l|377|. I3.12a|l 
and l|3.28|l . together with equation (|3.6|l imply that x(t) is of the form: 

x(t) = x + I 2 (t) +o(l). 



4 The presence of the term C e 2 2 7 in the bound of the second moment of the particle velocity 
can become important when studying the problem considered in this paper numerically. We refer to 
[Tgl for details. 
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Thus, in order to analyze the behavior of x(t) as e tends to we need to identify the 
contribution of the term I%{t) to the limiting equation. In this section we use the 
bounds on the moments of y(t) that we derived in section ?3. 31 to investigate precisely 

0. 



the limit of h(t) defined by (|3.5b(l as e 

We start with the regime 7 S (0, 2). We have the following: 

Lemma 4.1. Let x(t) be given by (13.1(1 and let y(t) = x(t) be given by 1(3. 2JI . 
Suppose that conditions 1(2. 9|l and ((2.12(1 with p — 1 hold. Then for 7 G (0, 2) the 
term I2 (t) given by (I3.5b(l has the form: 



whe 



hit) = / f(x(s))A- 1 dW(s) + Hit) 



E 



sup \\H{t)\\ 2 A < Ce {2 ~^ a , 

0<t<T J 



and where a > is arbitrarily small. 

Proof, hit) is given by equation ((3.24|) : 



(4.1) 



hit) 



fixis^A' 1 dWis) - eJx(t) + eJ 2 {t). 



We have that 



sup 

0<t<T 



iJ.it) 



2p 



< Ce 2 ^. 



Furthermore, estimate ((3.29(1 gives 



E [ sup 

\0<t<T 



eJ 2 (t) 



2p 



< £7 £ (2-7)p-er_ 



Estimate (14.1(1 follows from the above bounds. □ 

From our estimates on terms hit), i = 1, . . .4 we anticipate that xit) converges 
in mean square, as e — > 0, to Xit), which satisfies equation ((2.16(1 . The proof of this 
convergence is presented in section [3] 

Now we proceed with the case 7 £ [2, 00). We have the following lemma. 

Lemma 4.2. Let xit) be given by \ H. 1\) and let y(i) = ±(£) be given by 1(3.2(1 . 
Suppose that conditions 1(2.7(1 . 1(2.9(1 . ((2.10(1 and 1(2.13(1 hold. Then for 7 £ [2, 00) the 
term I2 it) in IS. 1\) has the form: 



hit) 



V • (fixi S ))Qf T ixis))) ds - / /(x(s))6V ■ fixis)) ds 



+ / fix(s))A~' dWis) 



c7"l 



df(x(s))y(s)Ti{s)ds + H(t), 



(4.2) 



whe 



sup \\H(t)\ 

0<t<T 



2p 



< Ct 2p - a . 
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where a > is arbitrarily small. 
Proof. From l|3.24|l we have: 

/"* °° f (x(s)) /X~ 
h{t)= 3 dp j {s)-eJ 1 {t)+eJ 2 {t), (4.3) 

Jo - =1 a j 

( - . . V 

with E sup < t <T 



eJi(t) J < Ce 2p a 5 . In order to study the term Ji(t) we need 
another two integrations by parts. We apply Ito formula to the function 

Gijk(x,y,i]) = ./',,> ;.r:///.n, \i, 

from which we obtain, after some algebra: 

fij,kVjVk dt = _ e i+i d ( f ',■ >■■ 'I i !!>■■ \ , e7 +i y- fjjMVjykVl dt 



£ 7+l^ /" fij,kVjyk \ f7 +i y^ fijMVj 

c fij,kVjbk dt + y^ fij,kfk P VjVp dt 



p=l a 3 



f . L:.l. !h\/\ d ^ _ e i-i f dt (4 .4) 

Now we define the following functions (no summation): 

Fijkp = fij,kfkp 

and 

Gijkp — FijkpT]p1]jQj 

We apply Ito formula to the function Gijkp to obtain: 

dGijkp = ^2 F ljkp ,£yeri]Vpaj 1 dt ~ -2 F ijkpaj 1 r l:j i]p(a p + a )dt 
1 e 

+ -^F ijkp aJ 1 y / X j \p5 : jpdt + -F i:jkp aj 1 {r] p ^/X' j dl3j + iy y^d/3 p ), 



from which, after multiplying through by e 2 (ctj + a p ) 1 and taking the sum over 
p = 1, 2, . . . we get: 

y^ fijjkfkpjhjh d ^ _ fijtkfkj^j j. 

/ 00 p j, \ d 00 / p p \ 

2^ / y^ jij^jkpVpVj \ £ 2 y^ y^ Uij,kjk P ),eye'np'nj dt 



-e 2 ' 



■DC 



-eJ2 T^T s bpV^Wi + mVKdPp)- (4-5) 
~i a iK a P + a o) 



5 This estimate is independent of 7. 
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Note that 



fij,kfkj {fijfkj)^ fijfk.kj- 

Thus, from the above calculations, after taking the sum over j <E Z d and k = 1, . . . , d 
we obtain: 

eJ 2 (i) = f V • (/(o;( S ))e/ T ( a; (5))) - /* /(x( S ))6V • / T (x( S )) d S 
Jo Jo 

+ [ f(x( S ))A- 1 dW(s)-e^- 1 f df(x(s))y( S )r,(s)ds + H(t), 
Jo Jo 

where H(t) — J2i=i He(t) with 



j,p=l k=l 



fij,k (x{s))f kp {x(s))r) P {s)Vj (s) 



(4.6a) 



m) = e* £ E 

j, P =ife,<=i 



(/i3,/c£(x(s))/fcp(x(s)))^y £ (s)?7p(s)?7 J (5) 
aj(aj+a p ) 



<is, (4.6b) 



r t oo d 

I ^ 

J ° j,p=lk=l 



fi3,k(x(s))fkp(x(s))v P ( S )^X] 
atjiotj+atp) 



d/3j{s), (4.6c) 



PL " 

= EE 

- 70 j>=ife=i 



' ' fijA x ( 8 ))fkp( x ( 8 ))Vj(.8)y/\> 



aj{aj + a p ) 



d/3 p (s), (4.6d) 



oo d 



fi 3 ,k(x(s))y k (s)r)j(s) 



3 = 1 k=l 



a,; 



(4.6e) 



H l 6 (t) = e~< +1 fy V fi3M{x{s))y k {s)yi(s)r) 3 {s) ^ 
Jo j= i M=1 a J 



(4.6f) 



„t oo d 

H M = * EE 

JO „-_i . 1 



fij,k{x(s))Vj( s )bk(x(s)) 



ds, 



j=i fc=i 



(4.6g) 



£ oo d 



#lw = ^/ EE 



fij,k( x ( s ))yk{s)y/% 



dfois). 



(4.6h) 
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Now we have to bound the terms Hg(t), £ = 1, . . .8. The necessary estimates are 
proved in Lemma I A . 1 1 which is presented in the Appendix. The final result is that 



E 



sup 

0<t<T 



\H(t)\ 



J 2p 



< Ce 2p ^. 



This completes the proof of the lemma. 



Lemma 14.21 together with the estimates on ^3 (t) , I4 (t) will enable us to show 
that for 7 e (2,oo) the particle position x(t) converges in L 2p (fl, C([0, T];R)) to the 
solution of the Ito SDE l|2.19|l . The precise convergence theorem will be proved in the 
next section. 

The above argument fails when 7 = 2 since in this case we cannot control the last 
integral on the right hand side of l|4.2|) uniformly in e. For the convergence theorem 
in this case we need the following corollary of the previous lemma: 

Corollary 4.3. Let the conditions of Lemma \4-^\ be satisfied and let 7 = 2. 
Then Li (t) has the following form: 



h{t) = 



V • 



(/(x( s ))6/ T (x( s ))) ds- ( f(x(s))QV-f T (x(s))ds 



f(x(s))A- 1 dW(s) + H(t), 



(4.7) 



with E sup 0<t<T H(t) 



2l> 



< e 2p a , where a > is arbitrarily small. 



Proof. We set 7 = 2 in 14. 4|) . solve for fij,kDk r lj OL j 1 dt and combine the result with 
(|4.5I) . sum over j S Z d and k — l,...,d and integrate over [0, t] to obtain 

h{t) = f n v ■ (/(x( s ))e/ T (x( s ))) d s - J f(x( s ))ev ■ f T (x( s )) ds 

f(x(s))A- 1 dW(s)+H(t), 

with H(t) — Sfe=i Hk(t). The terms Hk(t), k = 1, ... ,8 are similar to the terms 
defined in H4.6(l . with the difference that the diagonal operator (L + A)^ 1 is applied 

to every one of them. Similar techniques to the ones used in the proof of the previous 

2p^ 



ones still apply. We obtain estimate E ( sup 0<t<T 
that conditions l|2.13|l hold. □ 



H{t) 



< Ce 2p - a , provided 



The above corollary will enable us to show that, for 7 = 2, x(t) converges to 
X(t) which satisfies SDE (|2~2"T|) . This leads to the surprising conclusion that in this 
case the correction to the drift is not the usual Stratonovich correction. The precise 
convergence theorem will be proved in the next section. 

5. Proofs of The Convergence Theorems. In this section we prove the con- 
vergence Theorems 12.11 12.21 and 12.31 In the following proofs we will use the fact that 
/ : £2 —* R d is a Lipschitz continuous map, provided that condition (|2.12b|) with p — \ 
holds. In particular, a calculation similar to the one presented in the proof of Lemma 
13.41 equation (|3.15() yields: 



|| (f(X(s)) - f(x(s))) A- l \\ L , < C \\X(s) - x(s)\\ 



(5.1) 
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Proof of Theorem 12.11 We combine ijSJJ together with estimates (|3.7H , 
Corollary 13.101 and Lemma f4. II to write x(t) in the form: 

x{t)=x Q + ( f{x(s))A~ 1 dW(s)+ [ b(x(s)) ds + Ri, (5.2) 



with E (sup < t < T ||i?i|| 2p j < C ( e ( 2 -7)p-^ + e w- CT ). Now we take the difference be- 
tween X(t) given by H2.16fl and x(t) given by 1)5. 2|) . raise it to the 2p-th power, take 
the expectation value of the supremum, use the Burkholder-Davis-Gundy inequality 
and the Lipschitz continuity of f{x), b(x), together with the estimate on R\ to obtain: 

e( sup \\X(t)-x(t)f p ) <C [ T E\\(f(X(s)) - f(x(s)))A-X P ods 

\0<t<T J Jo 2 

r T 

+C 



sup \\X(t) - x(t)\\ 2p ) ds + CE[ sup ||i?i(i)|| 2p 

0<*<s / \0<t<T 

< C U 2 -~ ,)p -' 1 + e^ p - a ) + C [ E f sup \\X(t) - x(t)\\ 2 A ds. 

\ ' Jo \0<t<s J 



We apply now Gronwall's lemma to the above equation for the function 

|2p 



£(T) = E sup \\X(t) - x{ 

\0<t<T 

to conclude the proof of Theorem 12. II □ 

Now we proceed with the convergence theorems for 7 £ [2, 00). Let us consider 
the case 7 > 2. 

Proof of Theorem 12.21 . Consider the integral 

- ,7-1 



J(t)=e^ 1 / df(x(a))y{s)r,(s)ds. 
Jo 

We use Lemma f3 . 71 with p = to deduce 

E ( sup ||J(<)|| 2p ) < Ce 2p( -^ 2) - a . 

Now we combine l|3.1[) with estimates (|3.7H . <|3.26H . (I3.12a[) and Lemma f4. 21 to- 
gether with the above estimate, to write x(t) in the form: 



x(t)=x + B(x(s))ds+ f{x(s))A- 1 dW(s) + R 2 (t), (5.3) 
Jo Jo 

with E (sup < f < T ||i? 2 (i)|| 2 ) < C ( e 2p(7-2)-<T + e 2 P -^ and 

B(x(s)) = b(x(s)) + V • (f(x(s))<df T (x( S ))) - f(x(s))e\7 ■ f T (x(s)). 

Now, assumptions <|2.9|) and (|2.7|l imply that the drift term B(x) in equation 1)5. 3|) is 
Lipschitz continuous: 

\\B{X) -B(x)\\ < C\\X-x\\, (5.4) 
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under condition (|2.13c[l . The Lipschitz continuity of B(x), together (|5.1|l and the 
Burkholder-Davis-Gundy inequality give 



E ( sup - x(t)\\ 2 A < d ( e 2 P(-r-V-° + e 2p - a )+C 2 [ E f sup \\X(t) - x(t)\ 

\0<t<T / ^ > Jo \0<t<s 

from which Theorem 12 . 21 follows . upon applying Gronwall's lemma. □ 

Now we are ready to present the convergence proof and theorem for the case 
7 = 2. Since the proof is essentially the same as the one of Theorem 12. 21 we will be 
brief. 

Proof of Theorem[T3 We combine (pH|) with estimates IpTTfl . $FHfy . I|3.12a|) and 

Corollary 14. HI to write x(t) in the form-for 7 = 2: 

-l 



x(t)=x + B(x(s))ds + f{x(s))A- 1 dW{s) + R 3 {t), (5.5) 
Jo Jo 

with E (sup < t < T ||i? 3 (i)|| 2p ) < Ce 2p ~ a and 

B(x(s)) = b(x(s)) + V • (f(x( S ))Qf T (x( S ))) ~ f(x( S ))QVf T (x(s)). 



Assumptions assumptions H2.9JI . I|2.7|l and (|2. 13cfl and ensure that B(x) is Lipschitz 
continuous. As in the proof of the previous theorem, we take the difference between 
X(t) given by 1)2.21(1 and x(t) given by l|5.5|l ). raise it to the 2p-th power, take the 
expectation value of the supremum, use the Burkholder-Davis-Gundy inequality and 
the Lipschitz continuity of the terms in l|5.5l) to obtain. 

e( sup \\X(t)-x(t)\\ 2p ) <C 1 e 2p - a + C 2 [ T e( sup \\X {t) - x(t)\\ 2p ) ds, 

\0<t<T J Jo \0<t<T J 

Now we apply Gronwall's lemma to obtain estimate l(2.22|l . □ 
6. Applications. 

6.1. Inertial Particles in a Random Field. A model for the motion of inertial 
particles in turbulent flows was introduced in I^HlEnj. It consists of Stokes' law for 
the particle motion with the background divergence-free fluid velocity field being an 
infinite dimensional Ornstein-Uhlenbeck process. We assume that the motion takes 
place on the two-dimensional unit torus T 2 : 

tx = v(x, t) — x, (6.1a) 
« = V- L V, (6-lb) 

dijj = vAip dt + y^dW, (6.1c) 
where V -1 := (ir—,—7r—) T stands for the skew-gradient and ib denotes the stream 

v OX-2 ' OX\ ' Of 

function. Furthermore, W(x,t) denotes a Q- Wiener process on 



II := (/ e Ll er (T 2 )- f fdx = 
I ' Jr 
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Various asymptotic limits for (|6.1[1 were considered in |24| . Let us consider now the 
scaling limit considered in this paper. 

We assume that r = r e 7_1 , that the inverse noise correlation time v is of 0(e~ 1 ) 
and rescale time by t — ► t/e. Moreover, we expand the solutions of iltj.lcfl in terms 
of the eigenfunctions of the Laplacian on T 2 , e . We also set K = 2n2? \ {(0, 0)} 
and denote := {?/ G C ; r\k — V-k}, equipped with the standard £2 inner product. 
Setting To = 1 for notational simplicity, the rescaled equations (|6.1II can be written in 
the form: 

e ,^mim_ m (62a) 



f(x)£ = ik x e ik -*£ k , (6.2b) 

keK 



df] k = -\\k\ 2 Vk dt + -yfedfy, keK. (6.2c) 

with k x = [k 2 - ki] T , r\ = {r]k}kLi e & K ■ Moreover, {/?k(*)}*Li are mutually inde- 
pendent one dimensional standard Brownian motions satisfying the reality conditions 

Pk = 0-k- 

It was shown in [2H| that 

f®f = al with a =Y,^ (6-3) 

keK Z \ k \ 

where O is defined in equation (|2.18() and / stands for the identity matrix. A similar 
calculation reveals: 



fSf T = al with a = ^ 



; K 2\k\*(i + \k\2y 

Furthermore, the incompressibility of the velocity field implies that 

V • f = 0. 

The above calculations imply that, for the inertial particles problem whose motion is 
modelled by (|6.1|l . the ltd and Stratonovich interpretations of the stochastic integral 
coincide and that the limiting equation of motion is 



X(t)=x + f fiX^A-Uwis), 
Jo 



for all values of 7 > 0. In fact, Theorems 12.11 12"2l and 12 . 3l together with the properties 
of the eigenfunctions of the Laplacian on T 2 yield 

PROPOSITION 6.1. Let x(t) be the solution of (|6.2a() and 7 G (0,oo). Assume 

that 



VX^ < 00, for 7 G (0, 2) 



keK 
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and that 

^ \/~Xk\k\ < oo, for 7 e [2,oo). 



keK 



Assume further that conditions (|2.14() l|2.15|l hold and that the initial conditions for 
(I6.2c|l are stationary. Then x(t) converges, as e — > 0, to X{t) which satisfies 



t 

-1 



X(t) = x Q + / f(X(s))A- L dW(s), 







with 



E( sup \\X(t)-x(t)\\ 2p ) < C f e w + e( 2 -^ CT ) for 76 (0,2), 

\0<t<T / \ ' 



E( sup \\X(t)-x{t)\\ 2p ) < Ce 2 - a for 7 = 2 

\0<t<T 



and 



e( sup \\X{t)-x(t)\\ 2 A < C (e 2p -< 7 + e 2p (~<- 2 '>- a ) for 7 €(2, 00), 

\0<t<T / V ' 

where a > is arbitrarily small. The constant C depends on the moments of the 
initial conditions, the spectrum of the Wiener process, the operator A, the exponent 
p, the maximum time T and a. Physically we are looking at inertial particles in 
rapidly decorrelating velocity fields, over long times. The parameter 7 effects the non- 
dimensional mass of the particle, which is of C(e 7_1 ). Provided that 7 > 0, which 
includes a range of massive as well as light particles, the limiting particle motion 
is equivalent in law to a Brownian motion, see eq. (|(j.3fl . For 7 = 0, however, the 
limiting motion is that of the integrated OU process: the particle velocity is of OU 
type 

6.2. Diffusion in Solids. Consider now the motion of a particle in one dimen- 
sion under the influence of a periodic potential V(x), subject to dissipation: 

tx = -V'(x)-x. (6.4) 

We assume that the derivative of the potential can be written in the following Fourier 
sine series: 

00 

3=1 

We assume further that the control parameters \Xj are noisy and of the form 

(ij = /i" + -r]j(t/e 2 ), 



where are constants and are one-dimensional OU processes driven 

by mutually independent noises: 



drjj = —j rjjdt + \J \jd/3j, j = 1 ... 00 
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Substituting the above into (|6.4[) and assuming that the particle relaxation time r is 
of 0(e 7 ) we obtain 



_^ oo 

y x = — Vq(x) — x H — sin(ja;)r;j (t) 



(6.5a) 



3=1 

1 2 1 — 

dr/j = — ^ i Vjdt H — ^Xjdf3j, j = 1 . . . oo, 

where = — Sj!Li sin(ja;)^. We use now Theorems 12.11 l2~2l and l2~3l to deduce 

the following result. 

PROPOSITION 6.2. Let x(t) be the solution of (|6.5() and 7 S (0, 00). Assume that 



J2V^r'<^, for 7 6(0,2) 



and £/ia£ 



^i/X^oo, /or 7 € [2,oo). 
J'=l 

Assume further that conditions (|2.14(l (|2.15|) /ioW ; i/iai Vo(x) G C^(M) and £/iai the 
initial conditions for (|6.2cH are stationary. Then x(t) converges, as e — > 0, £0 X(i) 
which satisfies 



X(t) = { 



J 3 = L 3' 



with 



^ - vtf (*) + i Jo Eili Fife) sin ( 2 ^( s )) rfs + /J ££i ^ x ^) d &( s ) 

xo - V£(x) + \ J* ET=i F sin ( 2 ^( s )) ds + /o 1 sin(jX( S )) d/J,-(«) 



: 7 < 2 
: 7 = 2 
: 7 > 2 



E( sup || A"(i) — x(t)|| 2pN ) < C (e 7P + e^ 2-7 ^ -0 ") /or 7 

\0<*<T / \ ' 



e (0,2), 



E f sup ||X(t) - x(t)|| 2p ) < Ce 2 ^" / 

\0<t<T J 



or 7 



Ef sup ||X(t)-a;(t)|| 2p ) <c(e 2p - CT + e 2p ( 7 - 2 )- CT ) for 7 € [2, 00), 

\0<t<T / \ ' 

where a > is arbitrarily small. The constant C depends on the moments of the 
initial conditions, the spectrum of the Wiener process, the operator A, the exponent 
p, the maximum time T and a. We remark that for 7 e (0, 2) the particle motion is 
in the mean potential 

V (x) =E(V(x,t)). 

On the other hand, for 7 > 2, the limiting motion particle motion is in modified 
potential which depends discontinuously on 7 as 7 — > 2 + . 
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7. Conclusions. The Ito versus Stratonovich problem is studied in this paper 
for a class of infinite dimensional mean zero Gaussian random helds. It is shown that 
the correct interpretation of the stochastic integral in the limiting equation depends 
on the rate with which the particle relaxation time t p tends to 0, relative to that of 
the noise correlation time t„. In particular, it was shown that in the case where t p 
and r n tend to zero at the same rate, then the limiting stochastic integral in neither 
of Ito nor of Stratonovich type. 

The proof of our convergence theorems is based entirely on the pathwise tech- 
niques developed in j$] and used previously in |24) . rather than the weak convergence 
methods of e.g. Our techniques enable us to obtain strong, i.e. pathwise, conver- 
gence results as well as sharp upper bounds on the convergence rates. A drawback 
of the method employed in this paper is that it is applicable only for noise processes 
which can be expressed as solutions of stochastic differential equations, like the one 
used in this paper. In order to apply the results reported in this paper to more general 
classes of colored approximations to white noise, weak convergence techniques will be 
more appropriate. 

Appendix A. Estimates on terms Hi(t), i = 1,...8.. In this appendix we 
prove the following lemma. 

Lemma A.l. Consider the terms Hi(t), i = 1,...8 defined in (|4.6|) and set 
H(t) = Yn =i H i{t )- Assume that conditions (E7J), (|2.13a|) (EHJ), (EHHJ , <|2.12a|) with 
p = \ and (I2.12b|) with p = hold. Then the following estimate holds: 

E ( sup \\H(t)\\ 2p ) < Ce 2p . (A.l) 

\0<t<T / 



Proof. We start with H\(t). First we compute: 



EE 

j,p=l k=l 



fij,k (x(s))f kp (x(s) )rj p {s)r]j (s) 



p=l 



\p=i 
Ji{t) + j 2 (t). 



Now, calculations similar to the ones employed in the proof of Lemma 13.21 enable us 
to obtain 



E( sup || Ji(t)\\ 2p ) <GV 

\0<t<T 



sup ||J 2 WH 2P ) <Ce-°, 

0<t<T 



for a > is arbitrarily small, provided that conditions 12.12b|) with p = 1 and 12.12af> 
with p = 1 hold. We use the above estimate and the definition of Hi(t) to conclude: 



sup \\H 1 (t)f P ) <Ce^-\ 

0<t<T 



for a > is arbitrarily small. 



32 



G.A. PAVLIOTIS AND A.M. STUART 



We proceed now with H 2 (t). We define: 

j 

{f i j,M( x ( s ))fkp(x(s))) 4 y i (s)'q p (s)'q j {s) 



At) ■= E E 

j,p=i k,t=i 



aj(aj + a p ) 



Now we compute 

oo d 

nt) = E E 

j,p=i k,t=i 

oo ^S+r 

<^E^ 

3=1 



fij.kefX x ( s ))fk P (x(s))ye{s)r] p (s)ri :j (s) 
a j( a j + a p) 



E E 

j,p=l k,l=l 



fij,ke( x ( s ))fkp,£(x(s))ye{s)f] p (s)r] :j (s) 
a j ( a J + a p) 



~ a« +r Ms)\ , ^ a] +r \\y( S )\\\ V] ( S )\ ~ ^|^(«)| 



E 



^E 



cv. 



E~ p 



<c^|| y ( s )|| 2 



oo \ ^ / oo \ ^ 

E< r -%( S )ij + ^E a P +r_1 i^( s )ij 



for ( € 1. We use now calculations similar to those used in order to prove Lemma 
13.21 and estimate (I3.13|) . together with Lemma f3. 61 to deduce 



E sup \\H 2 {t)\\ 2p < C(e 2 ^- ff + e 4 P- 2 f?- ff ), 

\0<t<T J 

provided that conditions 12.12(1 with p = 1 hold. We choose now £ = 1 to obtain 

e( sup ||ff 2 || 2p ) <Ce 2p - CT . 

\0<t<T J 

Consider now the term Hiif). We introduce the cylindrical Wiener process 



P =i 



Now we can write H±(t) in the form 



where the map F 



H 4 (t) = e / FdW, 
Jo 



is defined as 



oo a 



/y,fc (x(s))f kp (x(s) )r)j (s) 



oy(ay +a p ) 



7 P , i = l,...d, V7 e ^ 2 
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We need to estimate the Hilbert-Schmidt norm of F. We have: 

oo d 



\ F \\ 2 L 2 (t 2 .«. d ) 



EE 

p=l i=l 



EE 

k=l j=l 



fijA x ( s ))fk P (x{s))r]j{s)^/\ p 



xl / 00 r, l3+r rv a + r 

<<?E E j p 



a ji a j + a p) 
2 



«5 ' a ph\VK 



< c E( a p +r_1 
P =i 



E< 

^=1 



/3+r-l 



< 



c E< 



/3+r-l 



%1 



provided that condition (|2.12a|) with p = ~ holds. We use now the Burkholder-Davis- 
Gundy and Holder inequalities, together with a calculation similar to the one used in 
the proof of Lemma T3. 21 to deduce: 



sup \\H 4 (t)\\ 2p )<Ce^ f E\\F\\ 2p ds<Ce 2p , 

<t<T J Jo V ; 



(A.2) 



provided that condition (|2.12b|) with p = 1 holds. Exactly the same analysis provides 
us with the estimate 



e( sup \\H s (t)fA < Ce 2 *>, 

\0<t<T / 



under conditions l|2.12b|l with p = i and (|2. 12a,|) with p = l. 
Now we consider term (t) . We have 

H 5 (t) = - e ^ 1 d/(x( S )) 2/ ( S )A- 1 , ? (5) 
Now, the calculations used in the proof of Lemma \'6. 71 together with l|3.13(l yield 



0<t<T 



E( sup \\H 5 (t)\\ 2p ) <Ce 2p ^ +1) E[ sup \\df (x(s))y(s)A- 1 T 1 {s) 



|2p 



0<t<T 



(A.3) 



for a > is arbitrarily small, provided that (|2.12b|l holds with p = 1. 

Now we proceed with H§(t). We use a simple variant of (|3.13Jl and assume 12.13a|l 
to compute: 

2p 



e( sup \\H 6 (t)\\ 2p ) <Ce 2p ^ [ T EY 

\0<t<T J Jo 



00 d 

EE 

j=l k,t=X 



fij,u( x ( s ))yk{s)yi{s)r)j{s) 



ds 



< (7 e 2p( 7 + 

< Ce 2p ^^ 1 ^ cr 
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where a > 0. 

Consider now H-?(t). This term can be written in the following form: 

H r (t)=e[ df(x(s))b(s)A~ 1 T 1 (s)ds. 
Jo 

We use the Holder inequality, together with a calculation similar to the one presented 
in the proof of Lemma T3. 71 as well as assumptions (|2.9I) . (I2.12a[l with p = 1 to obtain: 

e( sup \\H 7 {t)\\ 2 A <Ce 2p [ T E\\df(x{s))b{x{s))A- 1 7 ] (s)\\ 2p ds 

\0<t<T J Jo 

<Ce 2p . (A.4) 
Consider finally H%(t). We write it in the form 

H 8 (t) = e~< f df{x{s))y{s)A- 1 dW{s). 



We use the Burkholder-Davis-Gundy inequality to obtain 

sup \\H s (t)\\ 2 A <C^( [ T \\df(x(s))y(s)A- 1 f L0 ) . 

0<t<T J \Jo J 

Now we have 

\\df(x(s))y(s)A- 1 \\ 2 L o<C\\y(s)\\ 2 , 
provided that condition (|2.12b|) with p = holds. Thus: 

e( sup \\H 8 (t)\\ 2p ) <Ce 2 '< p E\\y{t)\\ 2p <Ce 2p{ ~<-V-< J . 

Putting now all the above estimates together we obtain QA.ip . □ 
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